Shot Noise in Magnetic Tunnel Junctions: Evidence for Sequential Tunneling 
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We report the experimental observation of sub-Poissonian shot noise in single magnetic tunnel 
junctions, indicating the importance of tunneling via impurity levels inside the tunnel barrier. For 
junctions with weak zero-bias anomaly in conductance, the Fano factor (normalized shot noise) de- 
pends on the magnetic configuration being enhanced for antiparallel alignment of the ferromagnetic 
electrodes. We propose a model of sequential tunneling through nonmagnetic and paramagnetic 
impurity levels inside the tunnel barrier to qualitatively explain the observations. 
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The discovery of the giant magnetoresistance |Tj fol- 
lowed by the observation of a large tunneling magnetore- 
sistance (TMR) at room temperature in magnetic tun- 
nel junctions (MTJ's) _3] have boosted interest in spin- 
dependent electron transport in magnetic nanostruc- 
tures, especially in the spin-dependent tunneling phe- 
nomena [1, During recent years, there has been a 
growing interest in controlling the TMR and also the 
statistics of tunneling events in MTJ's by nanostructur- 
ing of the insulating barrier Q . A variety of new electron- 
correlation mechanisms have been proposed, typically 
based on transport through double MTJ's with either 
an open or Coulomb-blockaded quantum dot (QD) con- 
tacted by ferromagnetic electrodes. Electric shot noise 
(SN) is a powerful tool for studying correlations of tun- 
neling processes in nanostructures beyond the capabili- 
ties of dc measurements @. The growing list of theoret- 
ically investigated topics regarding spin-dependent shot 
noise includes the noise asymmetry between parallel (P 
and antiparallel (AP) ferromagnetic (FM) alignment 
as well as continuous variation of the SN over the rela- 
tive angle between FM electrodes [8j], SN through an ar- 
tificial (QD) Kondo impurity contacted by magnetic 
leads, shot-noise enhancement by dynamic spin block- 
ade in tunneling through a small QD JJ)] , and shot noise 
for spin-polarized and entangled electrons with spin-orbit 
interaction in the leads [ll| . The scope of experimental 
efforts 0, has however so far been much more 
limited and inconclusive with regard to the nature of tun- 
neling electron correlations even in the conceptually sim- 
plest spintronic devices, viz. MTJ's, as manifested by 
shot-noise measurements. 

Current fluctuations due to discreteness of electron 
charge flowing through the structure out of equilibrium, 
which provide the shot noise, contain information not 
accessible by time-independent conductance. Sensitiv- 
ity to quantum statistics, interference, and interactions 
between electrons passing through the device has made 
SN an effective tool for investigating quantum trans- 



port in meso- and nano-structures [y]. In the absence 
of any correlations, Poissonian shot noise is practically 
frequency independent at low frequencies with the noise 
power given by S* = 2e/, in terms of the average current 
/. The Fano factor F = S/2eI representing normalized 
shot noise is in general lowered below 1 for noninteracting 
electrons due to fermionic statistics. Electron-electron 
interactions can either further suppress or enhance the 
Fano factor (even beyond the Poissonian value) . 

Despite the theoretical excitement about perspectives 
of using the shot noise for investigation of spin-polarized 
electrons, behavior of the SN even in simple nonstruc- 
tured MTJ's remains unclear. Jiang et al. [l^ re- 
ported an observation of the "full" SN (i.e., F ~ 1) 
in MTJ's with AP alignment of electrodes. Later the 
same group measured a strong suppression (down 
to F w 0.45) of the SN in magnetic tunnel junctions, 
which was not understood. Our Letter reports the first 
systematic investigation of the tunneling statistics in a 
magnetic tunneling device by measuring shot noise in 
Co(80 A)|Al203(14 A)|Py(100 A) MTJ's with and with- 
out Cr doping of the insulating barrier. We demonstrate 
a decrease of the Fano factor and its dependence on the 
alignment of the ferromagnetic electrodes for certain bar- 
rier conditions. 

Details of sample preparation have been published pre- 
viously [l^. For Cr-doped samples, the tunnel barri- 
ers were deposited in two steps. After deposition of 
the underlying Co electrode, a first tunnel barrier was 
formed by deposition and subsequent oxidation of 7-9 A 
of Al. Subsequently, sub-monolayer amounts of Cr were 
deposited on the AI2O3 surface, followed by a second 
Al layer deposition (5-7 A) and oxidation, resulting in 
a "^-doped" Al203|Cr|Al203 tunnel barrier. The noise 
measurements use a setup described in Ref. [l6| . which 
employs the cross-correlation method. This technique 
removes uncorrelated noise from the amplifiers and the 
noise of the leads and takes into account nonlinearity 
of the dynamic resistance while converting the obtained 
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FIG. 1: (Color online) Typical dynamic resistance obtained 
in P state for the Cr-free (a) and Cr-doped (b) junctions with 
0.8 A 5-layer measured at 300, 77, and 2 K. (c) Reduction 
of the TMR with applied voltage for Cr-free and Cr-doped 
MTJ's at r = 2 K. (d) Dependence of the ZBA(%) = 100 x 
[R{0 mV) - R{WO mV)]/R{WO mV) determined for the P 
alignment on Cr (at 2 K). 



voltage noise into current noise. Out of 13 samples in- 
vestigated, the shot noise was measured for 11 MTJ's: 
5 without and 6 with (5-layer of Cr in the middle of the 
barrier, ranging between 0.2 and 1.2 A in thickness. 

Figure [T] shows typical electron transport characteris- 
tics of the studied MTJ's. The dynamic tunneling resis- 
tance vs bias V [Figs. [lja),(b)] measured at three tem- 
peratures for P alignment proves pinhole- free MTJ's 
For all MTJ's studied, an asymmetric parabolic conduc- 
tance background [17] plus a zero-bias anomaly (ZBA) 
below T 77 K, appeared in the resistance of the junc- 
tion (Rj) [Figs. [Ha),(b)]. Presently, there exists several 
possible explanations of the ZBA's in MTJ's [H, [H, 
which consider magnon- or phonon- assisted tunneling or 
two-step tunneling through impurities inside the tunnel 
barrier which are also coupled to some additional degrees 
of freedom. Simultaneous ZBA and SN measurements 
on our samples suggest the ZBA is provided by sequen- 
tial tunneling through impurities accompanied with spin 
flips. 

Doping of the barrier with Cr enhances the normalized 
ZBA, although this trend presents rather large dispersion 
[Fig. md)]. Conductivity and TMR are generally sup- 
pressed when Cr thickness is increased, but the relations 
between these parameter and the nominal Cr concentra- 
tion are not strictly monotonic. We have found, how- 
ever, that the changes in the TMR are correlated with 
those of the tunneling resistance (see below). This can 
be understood as follows: As the barrier width and the 
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FIG. 2: (Color online) (a) Typical voltage noise measured for 
a Cr-free MTJ at T = 2 K with the applied currents (from 
bottom to up) of 3.4, 5 and 6 nA. (b) Voltage dependence of 
F on bias, for the Cr-free (filled) and 0.4 A Cr-doped (open 
stars), also measured at T = 2 K. The error bars show stan- 
dard deviations. 



resistance increase, the relative role of two-step tunneling 
increases, which generally reduces the TMR. The TMR 
is also monotonically reduced with the applied voltage 
both for the Cr-free and Cr-doped MTJ'sJsee Fig. [He)], 
in accordance with the previous reports [18| . 

The measured low-frequency noise has a typical form 
for MTJ's, with the 1// noise dominating at / < 100 Hz 
and the "white" noise dominating at / > 100 Hz. Fig- 
ure m^a) shows a typical voltage noise for the frequency 
and bias range where the 1// noise does not affect the 
data and the applied bias {eV » fcsT) ensures that SN 
presents the dominant contribution to the total noise. 
Fig. [DJb) shows a typical dependence of F on bias. For 
most of the undoped MTJ's, the Fano factor was reduced 
below the Poissonian value (F < 1), while for the Cr- 
doped MTJ's F was always close to one. 

Figure [3ja) shows the TMR and the Fano factor for 
the P alignment as a function of the resistance by area 
product {R X A) at T = 2 K. The Fano factor was aver- 
aged over the range 40 — 120 mV where it is nearly bias 
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FIG. 3: (Color online) (a) Dependence of the TMR and F (for 
the P state) on the resistance area product {R x A). Solid 
horizontal line marks F = 1. Vertical dashed line separates 
the Cr-free and Cr-doped regions, (b) Dependence of the 
relative variation of F with alignment AF/Fp{%) = 100 x 
(Fap - Fp)/Fp on the relative strength of the ZBA. Solid 
symbols point the undoped samples. The lines are guides for 
the eye. 
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independent. For the undoped MTJ's in the range where 
TMR is only weakly reduced with the product R x A 
(< 10^ Mfi/im^), we observed a gradual suppression of 
the Fano factor down to ^ 0.65. Doping of the barrier 
with Cr further increases the tunneling resistance and re- 
stores the Poissonian SN {F ^ 1). The suppression of F 
in a certain tunneling resistance range is not accompa- 
nied by the appearance of random telegraph noise as in 
Ref. [lj|, reduced TMR [2^, or by metallic temperature 
dependence R{T), clearly ruhng out pin-holes/hot spots 
across the barrier. Figurc[3I^b) shows the normalized AP- 
P F asymmetry as a function of the normalized ZBA for 
the P alignment. Surprisingly, we find that F depends on 
the alignment of the electrodes with .Fap /Fp > 1 only in 
the MTJ's with a weak zero-bias anomaly and becomes 
nearly independent of the alignment above some thresh- 
old value of the ZBA. We stress that the observed Fano 
factor asymmetry reflects only alignment of the FM elec- 
trodes, but not orientation of the magnetic field. 

Previous studies of the shot noise in non magnetic TJ's 
with AI2O3 barrier have observed Poissonian value F ~ 1 
[2l| . It was reported, however, that for non-magnetic 
TJ's with Si02 barrier, tunneling through localized states 
within the barrier, could indeed account for the measured 
reduced F (2^ . In the following, we consider two simple 
models for sequential tunneling via an island inside the 
tunnel barrier (see Fig. 3]), which capture some qualita- 
tive aspects of our measurements. First, consider tun- 
neling through a normal region {Rm) inside the tunnel 
barrier [Fig. dfa)]. Neglecting charging effects, we can 
simply sum the contributions to the (averaged) current 
and noise for the two spin species. To this end, suppose 
R M is coupled asymmetrically to the left and right reser- 
voirs (i?L and Rr) with the respective spin-dependent 
conductances given by 

9L]=9/\/P and gi^^^^g^/p, (1) 
9Li = ag/ y/]3 and gj^ = agy^ (2) 
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FIG. 4: Two models: (a) a large normal region Rm such that 
electrons can be treated as noninteracting, coupled to biased 
emitter and collector reservoirs in P or AP configuration, (b) 
7?M is a single spin-polarized impurity level that can hold only 
one electron with spin at a (random) angle 9 with respect to 
the collector magnetization. 




FIG. 5: Fano factor as a function of the left-right asymmetry 
parameter (3 setting a = 1/5 for the two models sketched in 
Fig. |4l Solid lines are for the P and dashed for the AP mag- 
netic configurations. Insets show the dimensionless current i 
and noise s = Fi defined in the text. Note the logarithmic 
scale for (3: Assuming tunnel rates depend exponentially on 
the barrier thickness, this corresponds to a linear scale for the 
Rm position inside the tunnel junction (small (3 correspond- 
ing to the proximity to Rl and large f3 to Rr). 



/3 is a dimensionless left-right asymmetry parameter and 
a characterizes spin polarization. The charge current at 
the voltage bias V is given by / = igV, parametrized 
by a dimensionless current i that depends on a and /3 
only. Let us furthermore write the zero-frequency shot 
noise as S" = 2esgV. The Fano factor thus becomes F — 
S/{2eF) = s/i. 

Let us recall first that, in general, noninteracting spin- 
less electrons in double-barrier structures have the series 
conductance and F (for spinless electrons) 

G = 9l9rI{9l +9r), {gl + gl)l{gL + gnf , (3) 

which are valid not only for large semiclassical Rm but 
also for sequential tunneling through a small Rm de- 
scribed by master-equation approach, in which case g's 
become respective transition rates instead of the tunnel- 
barrier conductances 6]. Summing corresponding cur- 
rent and noise for the two spin channels in the P configu- 
ration (neglecting correlations between two spin species), 
one trivially obtains for the current and F 



(l + «)V^/(l + /3), Fp = (l + /32)/(i + /3)2. (4) 



In the AP configuration. 



-Fap 



a(H 
a2(l 



a){l+P)^/[{a + fi){l 
+ 2/3- 2/32 + 2/33 + f3*) 



(a + /3)2(l-f a/3)2 



a/3)]. 



(5) 
(6) 



where (a <-> (3) is the same as the first summand but 
with a and f3 interchanged. These results are plotted in 
Fig.EKa) for a = 1/5. Note that Fap-Fp > for /3 - 1, 
roughly corresponding to the center of the junction for 
the Rm location, which is the region contributing the 
largest current, see Fig.[5Ka). The Fano factor asymmetry 
is reversed closer to the junction interfaces where the 
tunneling is asymmetric. 
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Consider now hopping through a single level that can 
hold only one extra electron, see Fig. IHb). If there is a 
large exchange-energy splitting along certain direction 0, 
one could imagine a situation when only spins polarized 
along 9 are energetically allowed to tunnel through. Wc 
can then calculate the current and noise using Eqs. ^ 
where the rates gL and now depend on 9 and the 
relative magnetic orientation in the leads [23| : 

5L =5lt(1 + cos0)/2 + .9l|(1-cos0)/2, (7) 
gn = 5flT (1 ± cos 9)/2 + (1 T cos 9)/2 (8) 

for the P (AP) configuration. Assuming 9 is ran- 
dom, we average the current and noise: — 
(1/2) J_^d{cos9).... This results in simple expressions 
for the P MTJ: 

tP - {l + a)^/[2{l+f3)] , Fp = ^ (9) 

which are the same as just averaging over 9 = and 
IT. There is no simple analytic form for the current and 
noise in the AP case. We plot the results in Fig. \5lh). 
Notice that the AP-P asymmetry is significantly reduced 
in comparison to Fig.O^a). 

In undoped MTJ's, we measured typically Fap > 
and both are significantly suppressed below 1, apart from 
the thinnest tunnel barrier. Both of these findings are 
consistent with the results in Fig. O for tunneling pre- 
dominantly through impurities in the middle of the bar- 
rier. Since F for tunneling through uniformly-distributed 
point-like localized states is in general 3/4 (in the absence 
of hopping correlations between the two spin species) 
[131 , which in particular applies to both models in Fig. IH 
the AP-P asymmetry would require some structural pref- 



erence towards tunneling through the middle of the bar- 
rier. The Fano factor is reduced to ~ 3/4, as shown 
in Fig. IHI^a), as the tunnel barrier becomes wider and 
the role of the two-step tunneling processes become rel- 
atively more important. The tunneling resistance does 
not indicate variable-range hopping involving multi-step 
tunneling, which was observed for wider tunnel barriers 
[1^. Observation of Poissonian noise after Cr doping 
could be due to an offset in Cr deposited nominally in 
the center of the junction, which leads to systematically 
asymmetric hopping. Finally, the observed correlation in 
the AP-P Fano factor asymmetry and the ZBA [Fig.[3]^b)] 
suggest that an inelastic spin-flip mechanism in the bar- 
rier is responsible for concurrent reduction of the former 
and enhancement of the latter. 

In summary, first systematic shot noise measurements 
in magnetic tunnel junction show an evidence for sequen- 
tial tunneling mediated by defects. We demonstrate for 
the first time that electron tunneling statistics can be 
manipulated by an applied magnetic field due to their 
dependence on the relative orientation of ferromagnetic 
electrodes and also by deliberately doping the tunnel bar- 
rier with impurities. Control over the sequential tunnel- 
ing could find applications in optimizing signal-to-noise 
ratio in magnetoelectronic devices and provide a new tool 
for investigating spin-dependent transport of electrons in- 
jected by ferromagnetic electrodes. 
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